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, A REVIEW  OF  RECENT  RESULTS  ON  SPREAD  F THEORY 


I.  Introduction 

Spread  F,  as  exhibited  by  diffuse  echoes  on  ionograms,  was 
discovered  some  forty  years  ago  by  Booker  and  Wells  (1938) . Until 
recently,  there  has  been  only  much  statistical  data  concerning 
Spread  F.  However,  in  the  last  few  years  with  advances  in  radar 
backscatter  measurements,  in  situ  measurements,  and  theoretical 
and  numerical  simulation  techniques  a clearer  picture  of  the  funda- 
mental plasma  instability  mechanisms  causing  equatorial  Spread  F 
phenomena  has  been  evolving.  It  should  be  emphasized  that  this  paper 
will  deal  only  with  equatorial  Spread  F theory  and,  in  addition,  dwell 
only  on  those  theories  using  plasma  mechanisms  as  a basis.  For 
experimental  results,  the  works  of  Farley  et  al.  (1970),  Dyson  et  al. 
(1974) , Kelley  et  al.  (1976) , Woodman  and  La  Hoz  (1976) , Morse  et  al. 
(1977) , and  McClure  et  al.  (1977)  should  serve  as  good  references  for 
the  interested  reader. 

In  speaking  about  Spread  F one  should  at  least  show  its  basic 

manifestation  on  an  ionogram  (see  Fig.  1) . After  all,  the  terminology 

Spread  F emerged  from  the  results  of  ionosonde  traces  such  as  those 

exhibited  in  Fig.  1.  The  multiple  traces  are  caused  by  magnetic  field 

aligned  irregularities.  If  there  were  no  irregularities,  a single 

crisp  trace  rising  to  the  right  would  be  exhibited  on  the  ionogram. 

At  this  point  we  want  to  remember  the  basic  equatorial  geometry. 

Figure  2 exhibits  the  basic  equatorial  nighttime  ionospheric  F region 

geometry,  i.e.,  the  geometry  under  which  equatorial  Spread  F occurs. 
Note:  Manuscript  tubmitted  November  16,  1978. 
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N (y)  represents  the  background  electron  density  as  a function  of 
altitude  (y) . Gravity,  £,  points  down,  the  ambient  magnetic  field,  B, 
is  horizontal  (pointing  north)  and  k represents  a horizontal  perturba- 
tion (in  the  westward  direction) . The  maximum  in  the  electron  density 
profile  is  the  F peak.  The  underside  of  the  profile  steepens  at  night 
due  to  chemical  recombination  effects  and  electrodynamic  forces.  The 
E region  has  been  severely  reduced  by  chemical  recombination  and  plays 
a negligible  role.  To  a plasma  physicist  this  geometry  is  a classi- 
cal flute  mode  geometry  and  one  might  expect  this  equatorial  iono- 
spheric geometry  to  be  unstable  to  a variety  of  plasma  instabilities. 
However,  until  recently,  one  of  the  basic  difficulties  was  getting  the 
unstable  irregularities  to  the  topside  (i.e.,  above  the  F peak)  when 
they  were  initiated  on  the  bottomside.  One  must  remember  that  the 
experimental  evidence  has  exhibited  both  top  and  bottomside 
irregularities . 

At  this  point  a general  brief  review  will  be  given  of  equatorial 
Spread  F (ESF)  theories.  First,  we  will  discuss  the  linear  theories. 
Dungey  (1956)  was  the  first  to  suggest  that  ESF  was  initiated  on  the 
bottomside  by  a Rayleigh-Taylor  instability.  Dagg  (1957)  suggested 
that  the  ESF  phenomena  was  due  to  E to  F region  coupling,  i.e., 
irregularities  in  the  E region  coupled  up  to  the  F region.  In  1959, 
Martyn  (1959)  was  the  first  to  suggest  that  ESF  was  a manifestation  of 
the  E X B gradient  drift  instability.  Calvert  (1963)  proposed  that 
the  downward  motion  of  the  neutral  atmosphere  at  night  was  responsible 
for  ESF.  This  mechanism  is  essentially  equivalent  to  the  E X B 
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instability  because  of  the  relative  motion  between  ions  and  neutrals 


in  determining  the  instability.  All  of  these  previously  invoked 


linear  instability  mechanisms  could  only  explain  the  formation  of 


bottomside  irregularities.  The  collisional  Rayleigh-Taylor  instability 


with  field  line  averaging  was  also  proposed  (Balsley  et  al.,  1972 


Haerendel,  1974)  as  a linear  instability  mechanism.  By  averaging 


(integrating)  the  density  along  the  magnetic  field,  the  total  electron 


raised  in  altitude  with  respect  to  the  local  electron  density  peak 


This  would  allow  the  linear  mechanism  to  operate  to  slightly  higher 


altitudes  ('v  100  km  greater)  , but  still  would  not  explain  the 


existence  of  irregularities  above  this  "new  peak".  Hudson  and 


Kennel  (1975)  pointed  out  the  importance  of  the  collisional  drift 


mode  in  ESF  in  the  wavelength  regime  30m  - 100m.  This  mode  could  be 


excited  on  both  the  top  and  bottomside  but  still  would  not  explain  the 


longer  wavelengths.  In  their  paper,  finite  larmor  radius  (FLR) 


Rayleigh-Taylor  instability 


Several  nonlinear  theories  have  been  invoked  to  explain  the 


different  ESF  observations.  For  example,  Hudson  et  al.  (1973)  sug 


gested  that  the  very  smallest  scale  (^.10m)  irregularities  (e.g 


process.  In  this  prescription  a longer  wavelength  instability  sets 


up  the  driving  conditions  for  the  shorter  wavelengths  to  become 


unstable.  This  is  similar  in  spirit  to  the  successful  two  step 


theory  (Sudan  et  al.,  1973)  proposed  for  Type  II  equatorial  E region 


electrojet  irregularities.  Haerendel  (1974)  suggested  that  the  range 

of  wavelengths  (many  kilometers  down  to  centimeters)  exhibited  by  ESF 

phenomena  was  due  to  a multi-step  process.  This  scenario  is  as 

follows:  (i)  the  collisional  Rayleigh-Taylor  (R-T)  instability  with 

horizontal  wavevectors  is  driven  by  gravity  and  the  background,  zero 

order  electron  density  gradient  scale  length  on  the  bottomside;  then 

(ii)  the  E X B gradient  drift  instability  with  vertical  wavevectors 

arises  due  to  the  horizontal  density,  large  amplitude  variations  set 

up  by  the  collisional  R-T  instability;  then  (iii)  the  inertia 

(collisionless)  dominated  R-T  instability  arises;  and  finally  (iv) 

kinetic  drift  waves  grow  upon  these  irregularities  after  they  reach 

-2 

large  amplitude.  Chaturvedi  and  Kaw  (1976)  interpreted  the  k 

measured  power  spectrum  of  the  ESF  plasma  density  irregularities  in 

terms  of  a two  step  theory.  In  this  theory  longer  wavelength  R-T 

modes  couple  to  kinetic  collisional  drift  waves  in  such  a manner  that 

-2 

the  mode  coupling  results  in  the  observed  k spectrum. 

I 

A major  breakthrough  was  made  by  Scannapieco  and  Ossakow  (1976) 
who  performed  a nonlinear  numerical  simulation  of  the  collisional  R-T 
instability  for  ESF  geometry.  The  simulation  results  showed  that  the 
collisional  R-T  instability  generated  irregularities  and  bubbles 
(plasma  density  depletions)  on  the  bottomside  of  the  F region  which 
subsequently  rose  beyond  the  F peak  by  nonlinear  polarization  induced 
E^  X B forces.  This  was  the  first  theoretical  result  to  explain  how 
long  wavelength  irregularities  could  appear  on  both  the  bottomside  and 
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topside  of  the  F region.  The  bubble  phenomena  was  in  accord  with  the 

recent  observations  (Kelley  et  al.,  1976;  McClure  et  al.,  1977, 

Woodman  and  La  Ho 2,  1976)  of  plasma  density  depletions.  An  analytical 

nonlinear  mode-mode  coupling  theory  for  the  coherent  development  of 

the  collisional  R-T  instability  was  performed  by  Chaturvedi  and 

Ossakow  (1977) . This  theory  suggested  that  vertical  modes  would  be 

-2 

dominant  and  result  in  a k power  spectrum.  Hudson  (1978)  extended 
the  previous  results  to  the  collisionless  R-T  regime  and  reached 
similar  conclusions.  Analytical  models  for  the  rise  of  collisional 
and  collisionless  R-T  ESF  bubbles,  in  analogy  with  fluid  bubbles, 
was  presented  by  Ott  (1978) . At  the  same  time,  Ossakow  and  Chaturvedi 

(1978)  presented  analytical  models  for  the  rise  of  collisional  R-T  ESF 

1 

bubbles  within  the  context  of  the  electrical  analogy  with  barium 
clouds. 

Costa  and  Kelley  ( 1978a, b)  suggested  that  coherent  steepened 

-2 

structures  and  not  turbulences  would  give  a k power  spectrum.  More- 
over, these  sharp  gradients  could  cause  small  scale  sizes  (^  20m)  by 
collisionless  low  frequency  (much  less  than  the  ion  gyrofrequency , ft^) 
kinetic  drift  waves  via  a two  step  process.  Their  analysis  was  a 
linear  one  carried  out  on  a nonlinear  state,  i.e.,  one  achieves  the 
steepened  gradients  by  nonlinear  processes  and  then  one  performs 
linear  theory  on  this  state.  Kelley  and  Ott  (1978)  suggested  that  the 
ESF  bubbles,  in  the  collisionless  R-T  regime,  generate  a wake  with 
vortices.  They  then  applied  two  dimensional  fluid  turbulence  theory 
to  the  model.  This  resulted  in  the  development  of  turbulence  at 
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shorter  and  longer  wavelengths  than  the  bubble  size.  This  in  turn  led 
to  a prediction  of  k ^ for  the  power  spectrum  (which  does  not  appear 
to  be  in  agreement  with  existing  experimental  observations)  in  the 
range  I,g  1 < k < LQ  1/  where  Lg  is  the  stirring  (bubble)  size  and  L^ 
is  a dissipation  length  cutoff.  In  a continuation  of  the  numerical 
simulation  work,  Ossakow  et  al.  (1978)  showed  a more  rapid  ESF  develop- 
ment and  higher  bubble  rise  velocities  resulting  from  sharper  bottom- 
side  background  electron  density  gradients  and  higher  altitudes  of  the 
F peak.  In  Huba  et  al . (1978)  very  small  scale  (wavelengths  *\<  lm  and 

I 

36cm)  irregularities  are  reported.  A two  step  process,  utilizing  high 
frequency  ,(^(K)  kinetic  drift  cyclotron  or  lower  hybrid  drift 
instabilities,  is  invoked  to  explain  them.  Linear  theory  for  these 
instabilities,  in  the  ESF  environment,  was  performed  on  the  nonlinear 
state  with  encouraging  results. 

The  above  introductory  remarks  glaringly  point  out  that  much  work 
and  significant  progress  in  the  theoretical  area  of  equatorial  Spread 
F has  been  accomplished  in  the  past  few  years  (indeed  just  look  at 
the  number  of  publications  during  1978  alone) . Notwithstanding  the 
recent  successes,  much  work  still  needs  to  be  done.  Indeed,  the 
theoretical  and  numerical  simulation  efforts  in  ESF  are  continuing 
along  a hot  and  heavy  path.  Section  II  of  this  paper  presents  out- 
lines of  some  of  the  theoretical  efforts  briefly  mentioned  in  the 
preceding  paragraphs.  Given  the  length  limitations,  it  would  be 
exceedingly  difficult  to  outline  all  of  the  theoretical  works 
mentioned  or  even  to  present  all  of  the  details  of  a few  works. 
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Section  II,  hopefully,  will  wetten  the  reader's  appetite  to  read  the 
referenced  works.  Section  III  presents  a summary  concerning  ESF 
theory . 


II.  Theory 

In  this  section  we  present  some  representative  theoretical  and 
numerical  simulation  works  with  the  appropriate  references. 

a.  General . The  basic  plasma  fluid  equations  applicable  to  the 
equatorial  Spread  F ionosphere  are  as  follows: 


r—  + V • (n  V)=P-vTn 
9t  a —a  R a 


V x B 

-T  Vn  - en  (-V<}>  + — — ) = 0 

e e e c 


m.n.  (tt  + V.  • V)  V.  = -T.  Vn. 
1 1 3t  — 1 —l  11 


V.  x B 

+ en.(-V<}>  + ) + m.n.  g - m.n.  v.  V. 

l c li—ii  m— l 


V • J = 0 


J = ne  (V.  - V,) 
— —i  — e 


In  the  above  equations  the  subscript  a denotes  species  (e  is  electron, 

i is  ion),  n is  density,  V is  velocity,  P is  the  production,  v is  the 

R 

chemical  recombination  rate,  T is  temperature,  V is  the  gradient 
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operator,  the  ambient  magnetic  field  (taken  to  be  uniform) , 

e is  the  electronic  charge,  c is  the  speed  of  light,  m is  mass,  £ is 
gravity,  is  the  ion-neutral  collision  frequency,  J is  current,  and 
the  electrostatic  approximation  has  been  made  where  E = 

Equation  (1)  is  the  continuity  equation,  (2)  and  (3)  are  the  electron 
and  ion  momentum  equations,  respectively,  (4)  is  the  divergence  of 

I 

the  current  and  (5)  is  the  current  equation.  What  we  have  in  mind  is 
to  apply  the  set  of  equations  (1)  - (5)  to  the  two  dimensions  perpen- 
dicular to  at  the  geomagnetic  equator,  making  various  approxima- 
tions. 

Assuming  a harmonic  perturbation  dependence  of  the  form 

exp-i  (k  • x - ut) , where  x denotes  perpendicular  to  B (k  is  hori- 
— —x  — o — 

zontal)  and  linearizing  equations  (1)  - (5)  we  obtain  for  the  linear 
growth  rate 


-h 


+ (v.  - 4 g • 

in  in 


V] 


/2  - 


a)  r a)  + iy 
r 


which  reduces  to 


— . — , v . >>  4 g/i 

, v.  L R in 
in 


(q)  2 

L - VR  ' V in<<:  4 9/L 


in  the  collision  dominated  and  inertia  dominated  regimes,  for 
kL>>  1.  Equations  (7a)  and  (7b)  represent  the  R-T  instability, 
including  recombination  damping,  in  the  respective  regimes.  (It 
should  be  noted  that  this  only  represents  instability  on  the  bottom- 
side  of  the  F region  where  the  first  term  in  (7a)  and  (7b)  is 
positive. ) 

b.  2D  Computer  Simulation  Results.  In  this  section  we  will 
outline  the  two  dimensional  ( J.  B^)  computer  simulation  results 
(Scannapieco  and  Ossakow,  1976;  Ossakow  et  al . , 1978) . The  simula- 
tions follow  the  nonlinear  evolution  of  the  collisional  R-T 
instability;  consequently,  in  eqn.  (3)  inertial  terms  (i.e.,  the  left 
hand  side  of  the  eqn.)  are  neglected.  Furthermore,  one  takes  for  the 
F region  « 1 (f2^  = eBQ/mic)  and  the  quasi-neutrality  assumption 

is  made,  i.e.,  n a n.  a n.  Equations  (1)  - (5)  then  become 
e l 


r—  + 7 • (n  V ) = 
d t a —a 


* VT  * 2 * 


- v_  (n  -n  ) 
R a ao 


V • J=0  , J = ne(V.  -V) 
— — — l — e 


where  one  has  taken  Tg  = T^  = 0 for  simplicity  (see  Ossakow  et  al., 
1978) , B = B z and  the  subscript  o refers  to  equilibrium  quantities. 


Making  the  electrostatic  assumption 


E = - V<J> 


and  breaking  the  potential  into  an  equilibrium  and  a perturbed 


quantity, 


= *0  + $ 


eqns.  (8)  - (11)  become 


_ £-(V$  x z)V  • n = - vn (n-n„) 


R o' 


V • (v  n V$)  = — (g  X z)  • Vn 

in  p *■ 


where  eqns.  (14)  and  (15)  are  taken  to  be  two  dimensional  (iB  ). 

- o 

Linearizing  eqns.  (14)  and  (15),  taking  a horizontal  perturbation 
results,  as  in  eqn.  (7a)  but  in  a more  illustrative  form,  in  the 
linear  growth  rate 

_2_.^o 

Y v . n VR 
in  o 


This  clearly  shows  that  only  the  bottomside  of  the  F region  where 
VnQ  is  positive  can  be  linearly  unstable  (and  only  if  the  first  term 

> V* 

Equations  (14)  and  (15)  were  solved  numerically  using  a 

vertical  mesh  spacing  of  Ay=  2 km  and  total  y extent  of  200km,  and 

an  east-west,  horizontal  mesh  spacing  Ax  = 200m  and  a total  horizontal 

extent  of  8km.  Realistic  profiles  of  v,  and  v_  as  a function  of 

in  r 
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altitude  were  utilized.  The  system  was  initialized  with  a 


direction  with  a wavelength  ^km  and  the  evolution  in  time  of  (14) 


and  (15)  was  followed  for  different  background  electron  density 


profiles.  Figure  3 shows  the  results  for  a background  electron 


profile  with  an  F peak  at  354km  and  a minimum  bottomside 


background  electron  density  gradient  scale  length,  L '''10km.  In  this 


case  at  t = 4000  sec  a bubble  (plasma  density  depletion)  is  clearly 


forming  and  beginning  to  rise  in  the  central  portion  of  the  mesh 


(note  n i n + n.).  The  isodensity  contours  are  such  that  the  maximum 


absolute  value  of  the  enhancement  or  depletion  is  in  the  center  and 


outer  contours.  At  t = 4000  sec  the  maximum  depletion  within  the 


bubble  is  54%  and  the  maximum  enhancement  over  the  mesh  is  84%.  At 


t = 8000  sec  we  notice  that  the  bubble  has  reached  the  altitude  of 


the  F peak,  with  the  innermost  contour  of  the  rising  bubble  repre 


senting  a 41%  depletion.  At  t = 10  sec  the  main  bubble  is  clearly 


through  the  F peak  with  an  innermost  depletion  contour  of  41%.  How- 


ever, in  the  ionosphere  below  the  bubble  near  x = 0 there  is  a 71% 
depletion  contour,  similarly  in  the  wings  near  |x|  = 4km.  The  inner 


most  enhancement  contour,  at  this  time,  represents  a 236%  enhancement 


with  a maximum  inside  this  contour  of  294%.  Note  that  the  top  of  the 


main  bubble  is  at  an  altitude  of  375km  while  the  bottom  trail  of  the 


bubble  is  at  an  altitude  of  270km.  Between  t = 8000  a’  9 10  sec  the 


bubble  has  risen  ''-24km  which  represents  a rise  velocity  'v-12m/sec 


Also  note  that  the  bubble  is  'Vlkm  wide.  Figure  4 depicts  contours 
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of  constant  induced  potential  ^ at  t - 10  sec.  This  shows  that 
the  more  isolated  part  of  the  high  altitude  bubble  depicted  in  Fig.  3 
is  acted  on  by  an  induced  electric  field  which  points  from  west  to 
east  and  is  dipolar  in  nature.  This  causes  the  bubble  to  rise  with  a 
(-c/Bq) V 4>^  x z velocity.  However,  the  lower  portion  of  the  mesh  is 
acted  on  by  an  induced  electric  field  '•'hich  points  from  east  to  west. 
This  field  is  much  weaker  than  the  induced  field  acting  on  the 
isolated  portion  of  the  central  bubble.  The  lower  altitude  electric 
field  causes  the  enhancements  and  depletions  to  move  downward.  Thus, 
the  lower  altitude  portion  of  the  central  bubble  becomes  captured  by 
the  enhancement? . 

Figure  5 presents  for  comparison  a case  in  which  the  background 
electron  density  profile  (shape)  was  kept  the  same,  but  the  entire 
profile  moved  up  in  altitude  so  that  the  F peak  was  at  434km.  All 
other  parameters  are  the  same  as  in  Figure  3,  except  and  vR  are 
taken  for  the  altitude  range  332km  to  532km  (those  used  in  this 
simulation) . One  can  immediately  note  the  more  rapid  time  evolution 
of  the  Spread  F process  with  respect  to  that  presented  in  Fig.  3. 

At  t = 700  sec,  a rising  bubble  with  an  innermost  contour  of  79%  and 
a maximum  depletion  inside  this  contour  of  84%  were  noted.  At 
t = 1000  sec  the  bubble  has  reached  the  peak  and  at  this  time  the 
innermost  depletion  contour  is  85%.  The  long  trail  associated  with 
the  high  altitude  bubble  has  a 100km  extension  to  lower  altitudes. 

At  t = 1400  sec  the  top  of  the  main  bubble  is  at  an  altitude  'vSOOkm 
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and  has  a long  trail  connecting  to  an  altitude  of  357km.  There  is  a 
maximum  70%  depletion  within  the  innermost  contour  of  the  high 
altitude  bubble.  Between  t = 103  and  1400  sec  the  top  part  of  the 
bubble  rose  'x'65km  and  this  represents  a rise  velocity  'V'160m/sec. 
Potential  contour  results  for  this  simulation  show  similar  patterns  to 
those  exhibited  in  Fig.  4.  Naturally,  the  induced  electric  fields 
causing  the  bubble  to  rise  in  the  present  case  are  stronger.  This 
spread  in  bubble  rise  velocities  has  been  observed  by  AE  satellite 
data  (McClure  et  al.,  1977) . 

Other  numerical  simulations  in  this  series  have  been  performed 
and  the  paper  by  Ossakow  et  al.  (1978)  should  be  consulted  for  more 
details.  The  basic  conclusions  reached  from  these  simulations  are  as 
follows:  (i)  the  collisional  R-T  instability  causes  linear  growth  on 

the  bottomside  of  the  equatorial  Spread  F region;  (ii)  plasma  density 
depletions  (bubbles)  steepen  on  their  top  and  nonlinearly  rise  to  the 
topside  ionosphere,  beyond  the  F peak,  by  polarization  (induced) 

E X B forces;  and  (iii)  high  altitude  of  the  F peak,  small  bottomside 
background  electron  density  gradient  scale  lengths,  and  large  percent- 
age depletions  yield  large  vertical  bubble  rise  velocities,  with  the 
first  two  conditions  favoring  collisional  R-T  linear  growth  (instabil- 
ity) . In  addition,  large  spatial  bubbles  with  similar  rise  velocities 
to  those  presented  here,  but  with  almost  100%  depletions,  have  been 
produced  by  numerical  simulations  (Zalesak  et  al.,  1978) . In  these 
cases  the  horizontal  mesh  covered  200km  in  extent  (intermesh  spacing 
Ax  = 5km)  and  a long  wavelength  (''-75km)  initial  perturbation  was 


13 


used.  Bubbles  with  horizontal  dimensions  'vSOkm  resulted.  The  NRL 


has  also  added  neutral  wind  effects  to  these  simulations  and 


found  that  an  eastward  neutral  wind  results  in  a westward  motion  of 


the  bubbles  in  addition  to  its  rise.  This  is  also  in  agreement  with 


1977) 


observations  of  bubble  motion  (see  McClure  et  al 


c.  Analytical  2D  Coherent  Mode  Coupling  Results.  A two 


instability  was  investigated  by  Chaturvedi  and  Ossakow  (1977)  using 


analytical  means  and  considering  coherent  mode  coupling  as  the 


saturation  mechanism.  Equations  (14)  and  (15)  were  utilized  with 


etc.  This  yields  the  following  coupled  nonlinear 


term  on  the  RHS  represents  damping,  and  the  second  term  on  the  RHS  is 


the  nonlinear  term.  Comparing  the  nonlinear  term  in  (16)  with  the 


last  term  in  (17)  results  in 


Therefore,  eqn.  (17)  is  treated  linearly  and  the  nonlinearity  is 


retained  in  (16) . A perturbation  of  the  form 


— = A.  , sin  (k  v-u>t)  cosk  x + A_  sin  2k  x 
n 1,1  y x 2,0  x 


is  chosen.  In  this  analysis  x is  vertical  (altitude)  and  y horizontal 
(east-west).  This  is  the  way  it  appears  in  the  reference.  For  the 
convenience  of  the  reader,  in  referring  to  the  original  work,  we  have 
kept  the  coordinate  system  of  the  reference. 


The  coupled  mode  equations  for  the  amplitudes  A1  , and  A_ 

1 # 1 2,0 


become 


^ Yi  i i ^ oi A i _ A _ « 

o t if  1 1 ; x l / 1 z , u 


2,0  a 7 

~3t  Y2,0  A2,0  + 2 Al,l 


where  the  coupling  coefficient  a = k k * g/k^  v.  , and  w = u + iv. 

x y in  r ' 

It  should  be  noted  that  the  linear  growth  rate  y is  negative,  i.e., 

2,0 

the  purely  vertical  mode  A is  linearly  damped.  Also  y.  is  posi- 

Z , U 1,1 

tive  and  represents  linear  growth  of  the  mode  A . In  the  saturated 
steady  state  one  has  3A  _/3t  = 3A_  _/3t  = 0 and  from  (19)  and  (20) 

J.  , 1 2 , U 

this  results  in 


2 k L 
x 


J Vi£ 

V k 2 gLk  2 

x y x 


) ■(>  ■-) 
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For  typical  values  of  the  parameters,  0 >:>  A1  1 and  shows 

-1  -2 
A « k (i.e.,  the  power  spectrum  would  k ).  This  represents  a 
2. 1 0 X X 

coherent  nonlinear  evolution  where  the  linearly  damped  mode  A is 
generated  nonlinearly  by  the  linearly  unstable  mode  A by  a har- 

j-  9 -L 

monic  generation  (A  contains  k and  A_  - contains  2k  ) . More 

1 , 1 X 2.  $ U X 

detailed  information  can  be  found  in  Chaturvedi  and  Ossakow  (1977) . 

d.  Analytical  Models  for  ESF  Bubbles.  First  we  will  discuss 

collisional  R-T  bubbles  (Ossakow  and  Chaturvedi,  1978)  in  which  the 

electrical  analogy  with  plasma  density  enhancement  (e.g.,  barium 

clouds)  was  utilized  and  results  obtained  for  general  2D  (J§0) 

bubble  shapes.  Here  eqns.  (14)  and  (15)  are  utilized  with  the 

further  simplifying  assumptions  of  neglecting  recombination  chemistry 

(v  = 0)  and  the  explicit  altitude  dependence  of  the  ion-neutral 
R 

collision  frequency,  . The  following  set  of  equations  result. 


- Y*(V$  x z)  • 7n  = 0 
o 


V • (nV$)  = E*  • Vn 


fi.  m. 


E*  = E + 


v.  e 
in 


g x z 


(23) 


(24) 


(25) 


where  an  ambient  horizontal  electric  field  E has  been  included  to 

— o 

show  more  generality.  Equation  (24)  can  be  thought  of,  in  general, 
as  a potential  equation  for  a dielectric  immersed  in  the  applied 
electric  field,  E*.  The  plasma  density  depletion  (bubble)  is 
analogous  to  the  case  of  a cavity  immersed  in  a dielectric  with  a 
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uniform  electric  field  (E*) . In  general  (23)  and  (24)  do  not  admit  of 
a two  dimensional  f jBq) , steady  state  solution.  However,  for  the  case 
of  a constant  density  inside  the  depletion  and  a constant  density  out- 
side the  depletion,  two  dimensional  steady  state  solutions  can  be 
obtained. 

For  a constant  density  depletion  with  an  elliptical  shape  one  has 
n(x,y)  = nQ  - ^(x.y) 


°D 


■ [‘  - ir FU'*1] 

*"  H[l  - 'I'"  - 'S’"] 

(1.  x > 0 

H(x)  = j 

* 0,  x < 0 


(26) 


(27) 


where  H is  the  Heaviside  function  (note  the  geometry  is  taken  such 
that  x is  vertical  and  y horizontal) . For  this  elliptical  shape, 


the  solution  to  (24) , neglecting  E , is 

~o 


3y 


fii  miqf  a 5n/no  \ 

® Y b + a(l  ~ Wn Q) ) 


(28) 


Using  -( c/B ) x z this  further  results  in  the  nonlinear  vertical 
o 

bubble  rise  velocity,  V , given  by 

B 


\n\h  + 


<5n/n 


a(l  - 6n/n  ) 
o 


) 


(29) 


Limiting  cases  of  (29)  for  sheet,  cylindrical  and  slab  bubble 
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geometries  are  given  by 


9/v. 

in 


(6n/n  ) 
o 

(1  - 

6n/nQ)  \ b « a 

(30a) 

(6n/nQ) 

(2  - 

6n/n  ) 1,  b = a 
o 

(30b) 

0 

, b <<  a 

(30c) 

For  typical  values  of  \>in  as  a function  of  altitude.  Fig.  6 exhibits 

rise  velocities  given  by  (30b)  for  various  percentage  depletions. 

Table  1 shows  the  rise  velocities,  in  units  of  g/v.  , for  various 

in 

shapes.  (Note  that  the  linear  case  comes  from  linearizing  (23)  and 

(24)).  All  of  the  geometry  results  can  be  expressed  in  a concise 
formula 


V -_SLf(^ 

B v.  n 
in  o 


(31) 


where  f(5n/nQ)  is  an  increasing  function  of  the  percentage 
depletion  6n/nQ.  Basically,  the  results  predict  that  high  altitudes 
and/or  large  percentage  depletions  yield  high  vertical  rise  velocities 
for  the  bubbles  (in  agreement  with  experimental  observations) . 

Collisional  and  collisionless  (inertia  dominated)  two  dimensional 
cylindrical  R-T  bubbles  have  been  studied  analytically  by  Ott  (1978) . 
This  study  is  based  on  the  analogy  with  fluid  dynamic  flows  and  brings 
forth  some  of  the  work  done  on  two-dimensional  fluids.  This  study 
begins  with  the  basic  equations  (1)  - (3) , considers  two  dimensions 
(oBq) , sets  ng  « n^%  n,  and  makes  the  assumption  that 


Di  ” St  • 


V,  v 


in 


(31) 


To  lowest  order,  using  (31),  one  obtains  from  (3)  a lqwest  order  ion 


18 


Equations  (34)  - (36)  form  a complete  set  of  equations 

(o)  'V 

sufficient  to  determine  the  unknown  quantities  , n,  and  p.  In 

the  limit  of  \>in  -*•  o*  eqns.  (34)  - (36)  are  identical  to  those  of 
an  ideal  incompressible  fluid.  At  this  point  the  philosophy  taken 
is  that  there  is  much  to  be  learned  concerning  bubbles  in  the  ESF 
ionosphere  from  the  extensive  studies  of  bubbles  in  fluids.  One 
then  uses  a stream  function  £ such  that 


V.  = z x V F 


After  a series  of  manipulations  one  can  obtain  an  equation  for  the 
stream  function.  In  the  collisional  and  collisionless  case  one 
finds  that  for  certain  values  of  thq  bubble  rise  velocity  £ will 
satisfy  the  equation,  for  a cylindrical  shape  (see  Ott,  1978  for 
more  details) . 

The  results  of  this  study  by  Ott  (1978)  predicts  the  bubbles 
to  be  cylindrical  (circular  cap  at  top)  in  two  dimensions.  The 
bubble  rise  velocity  in  the  collision-dominated  regime,  for  a 100% 
depletion,  is  given  by 

VD  = - — (38) 

B v. 
in 

which  is  altitude  dependent.  This  result  is  the  same  as  that  pre- 
dicted by  (30b)  for  6 n/nQ  = 1,  i.e.,  a 100%  depletion.  In  the 
inertia  (collisionless)  dominated  regime,  for  a 100%  depletion,  the 
bubble  rise  velocity  is  given  by 


VB  = S (Rg) 
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Irregularities . The  work  in  this  section  represents  essentially 
multilinear  calculations  using  kinetic  theory  resulting  in  plasma 


kinetic  drift  wave  modes.  It  is  multilinear  because  it  depends  on 
say  a two-step  process  whereby  linear  theory  is  performed  on  the  non- 
linear state.  The  driving  density  gradients,  in  these  calculations, 
are  thought  to  arise  through  a primary  instability,  driven  by  the  zero 
order  background  ionospheric  equatorial  F region  electron  density 
gradient,  achieving  a large  amplitude  state.  The  zero  order  iono- 
spheric electron  density  gradient  is  of  larger  scale  size  than  the 
primary  instability  electron  density  gradient  scale  size,  which  would 
be  of  the  order  of  the  instability  wavelength.  Because  the  calcula- 
tions are  kinetic,  they  employ  particle  distribution  functions. 

Kinetic  drift  waves  have  been  investigated  for  laboratory  plasma 
fusion  conditions  for  over  twenty  years,  so  a well  developed  formalism 
exists. 

Before  proceeding  with  specific  calculations  let  us  present  some 
general  concepts  regarding  kinetic  drift  waves  which  will  be  useful 
for  both  types  of  calculations  presented  in  this  section.  The  basic 
geometry  is  such  that 


B = B z 
— — o 


n = n (x) 
o o 


21 


(similar  to  the  zero  order  background  equatorial  ionospheric 


geometry)  where  here  nQ(x)  arises  due  to  the  primary  instability.  The 


orthogonal  coordinate  system  is  completed  by  the  ion  and  electron 
diamagnetic  velocities  being  along  the  y axis, 


V = (V..  - V,  ) y 
— d di  de  1 


(41a) 


VJ.  = (V.  /2ft.)  (31n  n /3x) 

di  11  o 


(41b) 


V.  = - (V  /2ft  ) (31n  n /3x) 

de  e e o 


(41c) 


where  V.  = (2T.  /m.  ) and  the  larmor  radius  is  defined  by 

i,e  i,e  i,e 


r =*  V /ft  . The  linear  analysis  is  then  performed  with  this  in  mind. 
a a a 


First  the  low  frequency  (u)  « ft^  collisionless  drift  wave 


calculations  of  Costa  and  Kelley  (1978a, b) , for  kr^  'v  l (r^  is  the  ion 


gyroradius) , will  be  presented.  These  calculations  were  meant  to  pro- 
vide a basis  within  which  to  try  and  account  for  the  3 meter  radar 
backscatter  observations  of  Jicamarca  (see  Woodman  and  La  Hoz,  1976) . 

A linear  kinetic  dispersion  relation  is  derived  using  perturbations  of 
the  form  exp  i |^k  • x - wt  j , where  k = k^  z + k y and  w = + iy . 

uiDtions  are  made  such  that  w « ft.  . V.  <<lui/k..  ! « \ 

e 


Furthermore,  assumptions  are  made  such  that  u)  « ft.,  V.  < <| /k,,  |«  V 

11  1 £ 


and  w>v.  , v.  , v..,  which  are  the  ion-neutral,  ion-electron,  and 
in  le  11  ' 


ion-ion  collision  frequencies,  respectively.  This  analysis  depends  on 
having  a k ^ , a component  of  the  wavevector  along  the  ambient  geo- 

magnetic field.  Figure  7 shows  some  growth  rate  results  from  these 
calculations.  The  growth  rate  is  in  units  of  the  ion  thermal  velocity 


1 


> 


* 
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divided  by  the  electron  density  gradient  scale  length 
(L  = nQ  (dnQ/dx)  . Figure  8 depicts  some  measured  inverse  electron 
density  gradient  scale  lengths  during  ESF  and  maximum  growth  rates  as 
a function  of  these  inverse  scale  lengths.  The  basic  results  of  these 
calculations  show  maximum  growth  for  k r.  a 1.5,  i.e.,  X s:  20m  (for 

-L  1 j. 

typical  ESF  parameters)  with  growth  rates  ^ 1 sec  1.  For  more 
detailed  analysis  concerning  these  low  frequency  drift  waves  applied 
to  ESF  see  Costa  and  Kelley  (1978a, b) . 

Now  we  present  the  high  frequency  drift  wave  analysis  of  Huba 

a 

et  al.  (1978) . In  this  reference  radar  backscatter  observed 
irregularities  with  wavelengths  of  1 meter  and  36  cm  at  Kwajalein  dur- 
ing ESF  conditions  are  shown.  In  an  effort  to  explain  these  very  short 
wavelength  irregularities  high  frequency  (to  ^ fl^)  drift  waves  were 
analyzed  for  ESF  conditions.  These  waves  are  the  so-called  drift 
cyclotron  (DC)  or  lower  hybrid  drift  (LHD)  instabilities  with  maximum 


growth  rates  for  k rg  'v.  1 (and  k^  = 0)  . No  k ,j  is  required  for 
these  instabilities.  The  parameter  determining  which  instability 


operates  in  a collisional  plasma  is 


Cf  = (vii/Qi)  (kxri)2  (42) 

Utilizing  the  linear  dispersion  relation  for  high  frequency  drift 
waves  for  Cf  <<  1,  instability  results  for 

L/r.  < U/21)  (m./m  ) ^ (43) 

i i e 

where  l is  the  harmonic  number  (u^  % ) . For  the  0+  ESF  ionospheric 

plasma  this  requires  the  electron  density  gradient  scale  length 
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L < 340m,  which  is  satisfied  (see  Fig.  8) . Growth  rates  for 


these  instabilities  are  given  by 


Y ss  (n>e/mi)  l «i 


(44) 


2 7 

However,  the  condition  C,  <<  1 implies  that  (k  r.)  n <<  2 x 10  and 

r ± 1 

for  k r % 1 this  means  n « 103  cm  3,  which  is  quite  restrictive, 
j.  e 

Longer  wavelengths,  i.e.,  smaller  values  of  k r ^ would  raise  the 
density  restriction  somewhat,  but  still  be  restrictive  for  ESF 


For  C,  ^ 1,  the  lower  hybrid  drift  instability  is  operative  and 


there  is  no  threshhold  condition  on  L.  Basically,  the  collisions 


which  increase  C 


instability  to  operate.  The  real  and  imaginary  part  of  the  fre 


quency  for  the  LHD  instability  are  given  by 


n the  collisionless  limit  it  should  be  noted  that  the  DC  instability 


transforms  into  the  LHD  instability  for  high  enough  ion  diamagnetic 


velocities  such  that 


which  for  ESF  conditions  implies  that  L ^ 30m.  Figures  9 and  10  show 


some  typical  results  from  the  analysis  of  the  high  frequency  drift 


The  results  of  this  analysis  predicts  that  the  lower  hybrid 


drift  instability  is  dominant  for  most  typical  ESF  ionospheric  para- 
meters. Also,  maximum  growth  of  the  instability  occurs  for 

k r n,  1 (X  <v  21  cm)  , although  good  growth  rates  can  occur  for 
j.  e x 

X 'x,  1 m.  Finally,  from  this  instability,  large  growth  rates 


x 

(y  Jh)  resulting  in  growth  times,  t = y less  than  a second  can 
occur.  For  more  details  of  this  work  see  Huba  et  al.  (1978) . 


III.  Summary 

Although  much  progress  has  been  made  in  the  theoretical  efforts 
directed  toward  the  equatorial  Spread  F ionosphere,  especially  in  the 
past  three  years,  more  has  to  be  done.  Also  the  burden  cannot  be 
placed  on  the  theoretician  alone.  Correlative  measurements  have  to 
be  made  prior  to  and  during  ESF  conditions.  It  is  not  sufficient  to 
make  a single  measurement  with  one  instrument  and  then  expect  a 
complete  theoretical  description  of  ESF.  One  needs  to  know  the  state 
of  the  ionosphere  with  respect  to  driving  parameters  such  as  back- 
ground electron  density  profiles  (bottomside  electron  density 
gradient  scale  lengths  and  height  of  the  F peak) , d.c.  electric 
fields,  neutral  winds,  and  ionic  mass  composition  in  order  to  build  a 
predictive  model.  In  addition,  in  order  to  compare  results  from  the 
predictive  model,  measurements  of  the  in  situ  fluctuating  component 
of  the  electric  field  and  plasma  density  have  to  be  made,  as  well  as 
radar  backscatter  measurements  of  the  very  small  scale  irregularities 
( 10m)  and  ground  measurements  of  satellite  signal  propagation 
amplitude  and  phase  scintillations. 
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Some  of  the  remaining  theoretical  problems  could  be  listed  as 
follows.  (1)  What  are  the  effects  of  changing  initial  conditions, 
including  ion  inertia,  and  including  neutral  winds  in  the  numerical 
simulations?  (2)  An  analytical  description  of  many  bubbles  rising. 

(3)  How  do  bubbles  decay  and  what  role  does  diffusion,  etc.  play? 

(4)  How  does  a turbulent  development  occur?  (5)  What  are  the  effects 

of  other  regions  of  the  ionosphere  (e.g.,  E region)  on  ESF?  (6)  A 

more  complete  study  of  collisional  effects  on  drift  waves  is  needed. 

What  determines  when  the  small  scale  irregularities  should  occur  (a 

more  quantitative  description)?  (7)  A determination  of  the  nonlinear 

saturation  of  small  scale  irregularities  (instabilities)  is  needed. 

(8)  What  are  the  effects  of  ? Indeed  points  (6)  and  (7)  are  tied 
to  the  more  general  question  of  what  is  the  relation  between  the  very 

small  scale  (^10m)  irregularities  (which  radar  backscatter  observes) 

and  the  longer  wavelength  fluid  type  (e.g.,  R-T)  irregularities  (which 

are  primarily  responsible  for  Spread  F seen  on  ionograms)?  I am  sure 

that  some  more  questions  and  points  could  be  raised.  However,  the 

above  list  should  keep  the  theoreticians  busy  for  a reasonable  time. 
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gure  was  taken  from  Hudson  and  Kenne 
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Fig.  3 - Contour  plots  of  constant  nj/nQ  for  the  simulation  with  an  F 
peak  at  354km  at  t * 2000,  4000,  8000,  and  10,000  sec.  The  small  dashed 
contours  with  a plus  sign  inside  and  the  solid  contours  with  a minus  sign 
inside  indicate  enhancement  and  depletions  over  the  ambient  electron 
number  density.  The  large  dashed  curve  depicts  the  ambient  electron 
number  density  (values  on  upper  horizontal  axis),  nQ,  as  a function  of 
altitude.  The  vertical  y axis  represents  altitude,  the  lower  horizontal 
x axis  is  east-west  range,  and  the  ambient  magnetic  field  is  along  the 
z axis,  out  of  the  figure.  Taken  from  Ossakow  et  al.  (1978) 


Fig.  U - Contours  of  constant  induced  (polarization)  potential, . 
over  the  mesh  (see  Fig.  3)  at  t = 10,000  sec.  Plus  and  minus  denote 
positive  and  negative  values,  with  values  decreasing  in  magnitude  as 
one  goes  from  the  innermost  to  the  outermost  contours.  The  large 
dashed  curve  is  nQ.  Taken  from  Ossakow  et  al.  (1978) 
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Fig.  7 - Growth  rates  as  a function  of  b for  the  low  frequency  collisionleas  drift  waves 
(Pi  is  the  ion  gyroradius).  The  curve  Y <<  ur  is  obtained  from  a small  growth  rate  approx 
imatlon  expression.  Taken  from  Costa  and  Kelley  (1978b) 
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Fig.  8 - Observed  Inverse  gradient  scale  lengths  (1/n)  (dn/dx)  as  a function  of  altitude 
during  the  Natal  rocket  flight.  Also  indicated  by  the  straight  line  is  the  maximum  growth 
ratecf  low  frequency  drift  waves  as  a function  of  (1/n)  (dn/dx).  Note  that  p.  is  the  ion 
gyroradlus.  Taken  from  Costa  and  Kf.lley  (1978a)  1 


Fig.  10  - Maximum  growth  rate  as  a function  of  ion  diamagnetic 
drift  velocity  for  the  DC  instability  (solid  line)  and  the  Ii!D 
instability  (dashed  line).  The  other  parameters  are  the  same 
as  in  Fig.  9.  Note  that  the  DC  instability  goes  into  the  IHD 
instbility  for  Vjj/Vi  ?>  0.11.  Taken  from  Huba  et  al.  (1978) 
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COMMAhOER 

NAVAL  SEA  SYSTEMS  COmAfO 
WASHINGTON,  O.C.  20362 

OICY  ATTN  CAPT  R.  PITKIN 


COFtWOER 

naval  space  surveillance  system 
OAHLGREN,  VA  22448 

OICY  ATTN  CAPT  j.  H.  BuRTQN 


department  of  energy 


department  of  energy 
ALBUQUERQUE  OPERATIONS  OFFICE 
P.  0.  BOX  5400 
ALBUQUERQUE,  NM  87115 

OICY  ATTN  Da  CON  FOR  D, 


OFFICER- IN-CHARGE 

naval  surface  weapons  center 

WHITE  oak,  SILVER  SPRING,  MD  20910 
OICY  ATTN  COOE  F*i 


JEPARTl^NT  OF  ENERGY 
LIBRARY  ROOM  G-042 
WASHINGTON,  D.C . 20545 

OICY  ATTN  OOC  CON  f OR  A.  LABOHITZ 


DIRECTOR 

STRATEGIC  SYSTEMS  PROJECT  office 
OEP ARTIST  OF  THE  SAW 
WASHINGTON,  D.C.  20376 
OICY  ATTN  NSP-2141 
OICY  ATTN  NSV-2722  FRED  WIMBERLY 


EG&G,  INC. 

LOS  ALAMOS  DIVISION 

P.  0.  BOX  809 

los  alamos,  nm  85544 

OICY  ATTN  Da  CON  FOR  J.  BREEDLOVE 


naval  SPACE  system  activity 
P.  0.  BOX  92960 
WORLDWAY  POSTAL  CENTER 
LOS  ANGELES,  CALl c . 90009 

01CV  ATTN  A.  B.  HAZZAR-D 

headquarters 

ELECTRONIC  SYSTEMS  DlVl SION/DC 
DEPARTMENT  OF  THE  AIR  FORCE 
HANSCOM  AFB,  MA  01731 

OICY  ATTN  OCKC  MAJ  J.  C, 


UNIVERSITY  OF  CALIFORNIA 
LAWRENCE  LIVERMORE  LABORATORY 
P.  0.  BOX  808 
LIVERMORE,  CA  94550 

OICY  ATTN  Da  CON  FOR  TECH  INFO  DEPT 
OICY  ATTN  Da  CON  FOR  l-389  R.  OTT 
OICY  ATTN  Da  CON  FOR  L-31  R.  HAGER 
OICY  ATTN  Da  CON  FOR  L-46  F.  SEWARD 


COmANOER 

FOREIGN  TECHNOLOGY  DIVISION,  AFSC 
WR1GHT-PATTERSON  AFB,  OH  45433 
OICY  ATTN  NICD  LIBRARY 
OICY  ATTN  ETDP  3.  BALLARD 


LOS  ALAMOS  SCIENTIFIC  LABORATORY 

P.  0.  BOX  1663 

LOS  ALAMOS,  **  87545 

OICY  ATTN  Oa  CON  FOR  J.  WOLCOTT 
OICY  ATTN  Da  CON  FOR  R.  F.  TASChcK 
OICY  ATTN  Da  CON  FOR  E.  JONES 
OICY  ATTN  Da  CON  FOR  u.  MALIK 
OICY  ATTN  Da  CON  FOR  R.  JEFFRIES 
OICY  ATTN  Da  CON  FOR  J.  ZINN 
OICY  ATTN  Da  CON  FOR  P.  KEATON 
OICY  ATTN  Da  CON  FOR  D.  wESTERVElT 


C0MMAM>£R 

ROFC  AIR  DEVELOPMENT  center,  afsc 
GRIFF ISS  AFB,  NY  1 34m 1 

OICY  ATTN  oa  LI3RARY/TSLD 
OICY  ATTN  aSE  V.  COTNE 


SAMSO/S Z 

POST  OFFICE  BOX  92960 
WORLD* Ay  POSTAL  CENTER 
LOS  AnGElES,  CA  90009 
(v> ACE  D€cENSc  SYSTEMS 
OICY  'Ytn  SZj 


SAWIA  LABORATORIES 
P.  0.  BOX  5800 
ALBUQUERQUE,  W 8*115 

OICY  ATTN  Da  CON  FOR  J.  MARTIN 

OICY  ATTN  Da  CON  FOR  W.  BROWN 

OICY  ATTN  Da  CON  FOR  A.  TMORNdAOUGn 

OICY  ATTN  Da  CON  FOR  T.  wRIGmT 

OICY  ATTN  Da  CON  FOR  D.  DAhlGREN 

OICY  ATTN  Da  CON  FOR  3141 

OICY  ATTN  Da  CON  FOR  SPACE  »ROwECT  DIV 


SANOIA  LABORATORIES 
LIVERMORE  LABORATORY 
P.  0.  BOX  969 
LIVERMORE,  CA  94550 

01CY  ATTN  OOC  CON  FOR  d.  MURPHEY 
01CY  ATTN  OOC  CON  FOR  T.  COC* 


BOEING  COMPANY,  THE 
P.  0.  BOX  3707 
SEATTLE,  wA  98124 

01CY  ATTN  G.  KEISTER 
01CY  ATTN  0.  MURRAY 
01CY  ATTN  G.  HALL 
01CY  ATTN  J.  KENNEY 


OFFICE  OF  MILITARY  APPLICATION 
DEPARTMENT  OF  ENERGY 
RASnINGTCN,  D.C.  20545 

QICY  ATTN  DOC  CON  FOR  0.  GALE 


CALIFORNIA  AT  SAN  DIEGO,  UNIV  OF 
IPAPS,  B-0I9 
LA  JOLLA,  CA  92095 

01CY  ATTN  H£MIY  G.  BOCXER 


OTHER  GOVERN  NT 


BROWN  ENG  I BERING  COH>ANY,  INC. 
CUMMINGS  RESEARCH  PARK 
HUNTSVILLE,  Al  55807 

OICY  ATTN  ROMEO  A.  DEL1BERIS 


central  intelligence  AGENCY 
ATTN  RO/ SI,  RM  5G48,  HQ  BLOG 
WASHINGTON,  D.C.  20505 

01CY  ATTN  OSI/PSIO  RM  5F  19 


CHARLES  STARK  DRAPER  LABORATORY 
555  TECHNOLOGY  SQUARE 
CAMBRIDGE,  MA  02159 

01CY  ATTN  D.  8.  COX 
01CY  ATTN  J.  P.  GILMORE 


DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
WASHINGTON,  D.C.  20234 

CALL  CORRES:  ATTN  SEC  OFFICER  FOR) 
01CY  ATTN  R.  .MOORE 


CCmjTER  SCIENCES  CORPORATION 
6565  ARLINGTON  3LVD 
FALLS  CHURCn,  VA  22046 
QICY  ATTN  H.  BLANC 
01CY  ATTN  JOHN  SPOOR 
01CY  ATTN  C.  NAIL 


DEPARTMENT  OF  TRANSPORTATION 
OFFICE  OF  THE  SECRETARY 
TAD-44.1,  ROOM  10402-3 
400  7TH  STREET,  S.W. 
WASHINGTON,  D.C.  20590 
01CY  ATTN  R.  LEWIS 
01CY  ATTN  R.  DOHERTY 


COMSAT  LABORATORIES 
LINTHICUM  ROAD 
CLARKSBURG,  MD  20734 
01CY  ATTN  G.  HYDE 


INSTITUTE  FOR  TELECOM  SCIENCES 
NATIONAL  TELECOMMUNICATIONS  & INFO  ADMIN 
BOULOER,  CO  80303 

01CY  ATTN  A.  JEAN  (UNCLASS  ONLY) 
01CY  ATTN  W.  UTLAUT 
01CY  ATTN  D.  CROC  IE 
01CY  ATTN  L.  BERRY 


CORNELL  UNIVERSITY 

DEPARTMENT  OF  ELECTRICAL  ENGINEERING 

ITHACA,  NY  14850 

01CY  ATTN  0.  T.  FARLEY  JR 


ELECTROSPACE  SYSTEMS,  INC. 

BOX  1359 

RIOORDSON,  TX  75060 

OICY  ATTN  H.  LOGS  TON 

01CV  ATTN  SECURITY  (PAUL  PHILLIPS) 


NATIONAL  OCEANIC  t ATMOSPHERIC  AOMIN 
cNVIROMENTAL  RESEARCH  LABORATORIES 
DEPARTMENT  OF  COWERCE 
BOULDER,  CO  80302 

OICY  ATTN  R.  GRUBB 

OICY  ATTN  AERONOMY  LAB  G.  REID 


ESL  INC. 

495  JAVA  DRIVE 
SUf#<YVALE,  CA  94066 

OICY  ATTN  J.  ROBERTS 
OICY  ATTN  JAMES  MARSHALL 
OICY  ATTN  C.  W.  PRETTIE 


DEPARTMENT  OF  DEFENSE  CONTRACTORS 


FORD  AEROSPACE  & CCmJNl CAT  IONS  CORP 
3939  FABIAN  WAY 
PALO  ALTO,  CA  94303 

OICY  ATTN  J.  T.  MATTINGLEY 


AEROSPACE  CORPORATION 

P.  0.  BOX  92957 

LOS  ANGELES,  CA  90009 

OICY  ATTN  I.  GARFUNCEL 
OICY  ATTN  T.  SALMI 
OICY  ATTN  V.  JOSEPHSON 
OICY  ATTN  S.  BOWER 
OICY  ATTN  N.  STOCKwELL 
OICY  ATTN  D.  OLSEN 
OICY  ATTN  J.  CARTER 
OICY  ATTN  F.  MORSE 
OICY  ATTN  SMFA  FOR  PWW 


GENERAL  ELECTRIC  COff*ANY 
SPACE  DIVISION 
VALLEY  FORGE  SPACE  CENTER 
GOOOARD  BLVD  KING  OF  PRUSSIA 
P.  0.  BOX  8555 
PHILADELPHIA,  PA  19101 

OICY  ATTN  M.  H.  BORTNER 


GENERAL  ELECTRIC  CO*»ANY 
P.  0.  BOX  1122 
SYRACUSE,  NY  13201 

OICY  ATTN  F.  RE1BERT 


analytical  SYSTEM  ENGINEERING  CORP 
5 OLD  CONCORD  ROAD 
BURLINGTON,  .1A  01803 

OICY  ATTN  RADIO  SCIENCES 


general  electric  co*»any 

TE^O-CENTER  FOR  AOVANCED  STUDIES 
816  STATE  STREET  (P.O.  DRAWER  QQ) 
SANTA  BARBARA,  CA  95102 
OICY  ATTN  DAS I AC 
OICY  ATTN  DON  CHANDLER 
OICY  ATTN  TOM  BARRETT 
OICY  ATTN  TIM  STEPHANS 
OICY  ATTN  WARREN  S.  <WP 
OICY  ATTN  WILLIAM  MCNAMARA 
OICY  ATTN  B.  GA«ILL 
01CT  ATTN  MACK  STANTON 


BERKELEY  RESEARCH  ASSOCIATES, 
P.  0.  BOX  985 
BERKELEY,  CA  94701 

OICY  ATTN  J.  WORKMAN 


icP€RAL  ELECTRIC  TECH  SERVICES  CO. , InC. 
nMES 

COURT  STREET 
SYRACUSE,  NY  13201 

01CY  ATTN  G.  MlLL-MAN 

GENERAL  RESEARCH  CORPORATION 
SANTA  BARBARA  DIVISION 
P.  0.  BOX  6770 
SANTA  BARBARA,  CA  931 II 

01CY  ATTN  JOHN  ISE  JR 
01CT  ATTN  JOEL  GAR0AR I NO 

GEOPHYSICAL  INSTITUTE 
UNIVERSITY  OF  ALASKA 
FAIRBANKS,  AT.  99701 

CALL  CLASS  ATTN:  SECURITY  OFFICER) 
01CY  ATTN  T.  N.  DAVIS  (UNCL  ONLY) 

01CY  ATTN  NEAL  BROWN  (UNCL  ONLY) 

OICY  ATTN  TECHNICAL  LIBRARY 

GTE  SYL VANIA,  INC. 

ELECTRONICS  SYSTEMS  GRP-EASTERN  DIV 
77  A STREET 
NEEDHAM,  MA  021 94 

OICY  ATTN  MARSHAL  CROSS 

ILLINOIS,  UNIVERSITY  OF 

department  of  electrical  engineering 
JRBANA,  IL  61803 

OICY  ATTN  K.  YEH 

ILLINOIS,  UNIVERSITY  OF 
107  COBLE  HALL 
801  S.  WEIGHT  STREET 
URBANA,  IL  60680 

(all  corres  attn  security  supervisor  for) 

OICY  ATTN  K.  YEH 


<AMAN  SCIENCES  CORP 
P.  0.  BOX  7463 
COLORADO  SPRINGS,  CO  80933 
OICY  ATTN  T.  MEAGHER 

L INC ABI  T COR? 

10453  ROSELLE 
SAN  DIEGO,  CA  92121 

OICY  ATTN  IRWIN  JACOBS 

LOWELL  RSCH  FOUfCATION,  UNIVERSITY  OF 
450  AIKEN  STREET 
LOWELL,  MA  01854 

OICY  ATTN  K.  BIBL 

M.l.T.  LINCOLN  LABORATORY 
P.  0.  BOX  73 
LEXINGTON,  MA  021 73 

OICY  ATTN  DAVID  M.  TOWLE 
OICY  ATTN  P.  WALDRON 
OICY  ATTN  L.  LOUGHLlN 
OICY  ATTN  0.  CLARK 

MARTIN  .MARIETTA  CORP 
ORLANDO  DIVISION 
P.  0.  BOX  5837 
ORLAhOO,  FL  32805 

OICY  ATTN  R.  HEFFNER 

MCDOWELL  DOUGLAS  CORPORATION 
5301  30LSA  AVENLE 
HUNTINGTON  BEACH,  CA  92647 
OICY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
OICY  ATTN  GEORGE  MR02 
OICY  ATTN  W.  OLSON 
OICY  ATTN  R.  *.  HALPRIN 
OICY  ATTN  TECHNICAL  LIBRARY  SERVICES 


INSTITUTE  for  DEFENSE  ANALYSES 
400  ARMY -NAVY  DRIVE 
ARLINGTON,  VA  22202 

OICY  ATTN  J.  M.  AE IN 
OICY  ATTN  ERNEST  BAUER 
OICY  ATTN  MANS  WOLFMARD 
OICY  ATTN  JOEL  BENGSTON 

MSS,  INC. 

2 ALFRED  CIRCLE 
BEDFORD,  MA  01730 

OICY  ATTN  DONALD  HANSEN 

INTL  TEL  i TELEGRAPH  CORPORATION 
500  WASHINGTON  AVENUE 
NUTLEY,  NJ  07110 

OICY  ATTN  TECHNICAL  LIBRARY 

JAYCOR 

1401  CAMI NO  DEL  MAR 
DEL  MAR,  CA  92014 

OICY  ATTN  S.  R.  GOLDMAN 


JOWS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOWS  HOPKINS  ROAD 
LAUREL,  *>  20810 

OICY  ATTN  DOCJMENT  LIBRARIAN 
OICY  ATTN  THOMAS  POTEMRA 
OICY  ATTN  JOW  DASSOULAS 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 
OICY  ATTN  P.  FISCHER 
OICY  ATTN  W.  F.  CREVIER 
OICY  ATTN  STEVEN  L.  GUTSCHE 
OICY  ATTN  D.  SAPPED IELD 
OICY  ATTN  R.  BOGUSCH 
OICY  ATTN  R.  HENDRICK 
OICY  ATTN  RALPH  KlLfl 
OICY  ATTN  DAVE  SOWLE 
OICY  ATTN  F.  FAJEN 
OICY  ATTN  M.  SCME10E 
OICY  ATTN  CONRAD  L.  LONGMIRE 
OICY  ATTN  WARREN  A.  SCHLUETER 

MITRE  CORPORATION,  THE 
P.  O.  BOX  208 
BEDFORD,  MA  017  30 

OICY  ATTN  JOW  .MORGANS TERN 
OICY  ATTN  G.  HARDING 
OICY  ATTN  C.  E.  CALLAHAN 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 

OICY  ATTN  W.  MALL 
OICY  ATTN  W.  FOSTER 


LOCXtCED  MISSILES  t SPACE  CO  INC 
P.  0.  BOX  504 
SUWYVALE,  CA  94088 

OICY  ATTN  DEPT  60-12 
OICY  ATTN  0.  K.  CHURCHILL 

LOOCXED  MISSILES  A tO  SPACE  CO  INC 
3251  HANOVER  STREET 
PALO  ALTO,  CA  94304 

OICY  ATTN  MARTIN  WALT  DEPT  52-10 
OICY  ATTN  RICHARD  G.  JOWSON  DEPT  52-12 
OICY  ATTN  W.  L.  I*KDF  DEPT  52-12 
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PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVE RF I ELD  BLVD. 

SANTA  MONICA,  CA  90404 

OICY  ATTN  E.  C.  FIELD  JR 

PEWS YL VANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 

(NO  CLASSIFIED  TO  THIS  ADORESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 


l 


PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  .**  02173 

OiCY  ATTN  IRVING  L.  KOFSKY 

PHYSICAL  DYNAMICS  INC. 

P.  0.  SOX  J027 
BELLEVUE,  WA  98009 

01CY  ATTN  E.  J.  FREMOUW 

PHYSICAL  DYNAMICS  INC. 

P.  0.  BOX  1069 
BERKELEY,  CA  9**701 

01CY  ATTN  A.  TnOMPSON 

RIP  ASSOCIATES 

P.  0.  BOX  9693 

MARINA  DEL  REY,  CA  90291 

01CY  ATTN  FORREST  GILMORE 
01CY  ATTN  BRYAN  GA86AR0 
01CY  ATTN  WILLIAM  B.  WRIGHT  JR 
01CY  ATTN  ROBERT  F.  LELEVIER 
01CY  ATTN  WILLIAM  J.  KARZAS 
01CY  ATTN  H.  ORY 
01CY  ATTN  C . MACDONALD 
01CY  ATTN  R.  TuRCO 

RAM)  CORPORATION,  Th£ 

1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 

01CY  ATTN  CULLEN  CRAIN 
OiCY  ATTN  ED  BEDROZIAN 

RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

OiCY  ATTN  VINCE  TRAPANI 

SCIENCE  applications,  INC. 

P.  0.  BOX  2351 
LA  JOLLA,  LA  92038 

OICY  ATTN  LEWIS  M.  LlNSCN 
OICY  ATTN  OANIEL  A.  HAM. I N 
OICY  ATTN  D.  SACHS 
OICY  ATTN  E.  A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGALL 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

OICY  ATTN  BARBARA  ADAMS 

SCIENCE  APPLICATIONS,  INC. 

HUNTSVILLE  DIVISION 
2109  W.  CLINTON  AVENUE 
SUITE  700 

HUNTSVILLE,  AL  35805 

OiCY  ATTN  DALE  H.  DIvIS 


SCIENCE  APPLICATIONS,  INCORPORATED 
8400  WESTPARK  DRIVE 
MCLEAN,  VA  22101 

OICY  ATTN  J.  COOCAYf* 

SCIENCE  APPLICATIONS,  INC. 

80  MISSION  ORIVE 
PLEASANTON,  CA  94566 
OICY  ATTN  S2 


SRI  INTERNATIONAL 
333  RAVENS WOOD  AVENUE 
MEN.0  PARK,  CA  94025 

OiCY  ATTN  DONALD  NEILSON 
OICY  ATTN  ALAN  BURNS 
OICY  ATTN  G.  SMITH 
OICY  ATTN  L.  L.  COBB 
OICY  ATTN  DAVID  A.  JOHNSON 
OICY  ATTN  WALTER  G.  CHESNUT 
OICY  ATTN  CHARLES  L.  RINO 
OICY  ATTN  WALTER  JAYE 
OICY  ATTN  M.  BARON 
OICY  ATTN  RAY  L.  LEAOABRAND 
OICY  ATTN  G.  CARPENTER 
OICY  ATTN  G.  PRICE 
OICY  ATTN  J.  PETERSON 
OICY  ATTN  R.  HAKE,  JR. 

OICY  ATTN  V.  GONZALES 
OICY  ATTN  D.  MCDANIEL 


TECHNOLOGY  international  CORP 
75  WIGGINS  AVENJE 
BEDFORD,  F*  01730 

OICY  ATTN  W.  P.  BOQUIST 

TRW  DEFENSE  C SPACE  SYS  GROUP 

ONE  SPACE  PARK 

REDONDO  BEACH,  CA  90278 

OICY  ATTN  R.  K.  PLEBUCH 
OICY  ATTN  S.  ALTSCHULER 
OICY  ATTN  D.  DEE 

V1SIDYNE,  INC. 

19  THIRD  AVENUE 

NORTH  WEST  INDUSTRIAL  PARK 

BURLINGTON,  MA  01803 

OICY  ATTN  CHARLES  HUMPHREY 
OICY  ATTN  J.  W.  CARPENTER 
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